Fission of a multiphase membrane tube 
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A common mechanism for intracellular transport is the use of controlled deformations of the 
membrane to create spherical or tubular buds. While the basic physical properties of homogeneous 
membranes are relatively well-known, the effects of inhomogeneities within membranes are very much 
an active field of study. Membrane domains enriched in certain lipids in particular are attracting 
much attention, and in this Letter we investigate the effect of such domains on the shape and fate 
of membrane tubes. Recent experiments have demonstrated that forced lipid phase separation can 
trigger tube fission, and we demonstrate how this can be understood purely from the difference in 
elastic constants between the domains. Moreover, the proposed model predicts timescales for fission 
that agree well with experimental findings. 

PACS numbers: 87.16.Dg , 87.16.Ac , 68.03.Cd , 68.47.Pe 



Internal organization is one of the most intriguing as- 
pects of the cell. Living cells have to actively maintain 
gradients of all sorts. Compartment alisation and traf- 
ficking aid it in doing so, and both processes extensively 
use membranes. Not only are the various organelles in 
eukaryotic cells surrounded by membranes, but the ba- 
sic intermediates in the intracellular transport pathways 
as well are membrane structures such as tubes and vesi- 
cles 0. The generation and properties of these struc- 
tures have been extensively studied, and much is already 
known about their biology, biochemistry |2| and their 
biophysics 0, H • The emerging view is that the shape 
of the bilayer membrane in vivo is controlled not only by 
embedded and associated proteins p but also to a large 
extent by the mechanical properties of the bilayer itself 
0,0. For tubular structures in particular, mechanical 
effects play a major role: recent biomimetic experiments 
^ have shown that kinesin motors walking on micro- 
tubules can exert pulling forces on the membrane and 
prompt the formation of membrane tubes that resemble 
tubules identified in living cells. 

The existence of small membrane domains with a lipid 
composition that is markedly different from that of the 
rest of the membrane (sometimes referred to as "rafts" 
although considerable debate remains as to their precise 
interpretation) appears to be another key element of in- 
tracellular vesicular traffic Q , and also seems to be impli- 
cated in a multitude of cellular processes [8| . The hetero- 
geneity in membrane composition can be attributed to a 
phase transition leading to a local segregation between 
the various lipids constituting the membrane 9]. Sph- 
ingolipid domains in particular have been shown to be 
more structured than a classical liquid membrane due to 
specific interactions between their constituents 0. Un- 
der appropriate conditions they tend to aggregate into 
so-called liquid- ordered domains which are mechanically 
stiffer than the rest of the bilayer. Recently, an experi- 



FIG. 1: Breakage of a heterogeneous membrane tube [T3| . 
The brighter (and thinner) section at the tip on the left is 
a liquid-disordered DOPC domain. Fission events occur at 
the sites of formation of small domains resulting from phase 
separation. The time between two consecutive pictures is one 
second. Scale bar, 10 /xm. 



mental model system of vesicles including "raft-like do- 
mains" has been developed jlCj : it provides an elegant 
and efficient tool to study their properties in a more con- 
trolled way than in vivo. This procedure allows for sys- 
tematic studies of the effects of membrane composition 
[l^ . temperature changes 9] and protein absorption on 
the domain [~l| . 

The physics of membrane tube formation from ho- 
mogeneous vesicles has been studied both theoretically 
|l4l Ljjj and experimentally . Recent experiments in- 
volving one of us 16] study the interplay between lipid 
domains and the behavior of tubes, by pulling tubes from 
model membranes. Fig. illustrates one of the suprising 
conclusions of these experiments - a sequence of snap- 
shots taken at regular intervals (one second between two 
pictures) show an initially homogeneous tube that first 
undergoes phase separation (triggered experimentally by 
photoinduced oxidization of cholesterol), and, after about 
one second, ruptures precisely at the phase boundary 
and disconnects. The two lipid phases are easily dis- 
tinguished, once separation has occurred, by the use of a 
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fluorescent marker that preferentially sits in the liquid- 
disordered domains. Furthermore, the same experiments 
show that fission events such as these happen only in 
the phase separated tubes - tubes in which the lipids are 
mixed are essentially stable indefinitely. 

Statement of problem and summary. In this Letter, 
we address the dramatic loss of stability following phase 
separation from a mechanical point of view. We extend 
the theoretical models developed for homogeneous tubes 
[l^ to study the junction between two distinct phases, 
each of which far away from the junction has a tubu- 
lar shape. Experiments suggest that phase separation 
occurs on a much faster timescale than fission, and that 
the nucleation of the two phases leads to the formation of 
cylindrical domains between a more rigid and a less rigid 
phase. We therefore choose not to model the dynamic 
of the phase separation process [HI . The tube radii and 
the junction length are generally small compared to the 
length of each phase domain. In order to minimize the 
inter facial energy between adjacent domains, the inter- 
face rapidly becomes a circle perpendicular to the tube 
direction. The coarsening stage of the phase separation 
process proceeds very slowly to eventually form two ho- 
mogeneous phases in equilibrium, but this slow relaxation 
is always preempted by tube fission. 

We assume here that the tube and junctions are axi- 
symmetric with respect to the direction along which the 
tube is pulled (the z axis). We consider one junction be- 
tween two semi-infinite tubes each consisting exclusively 
of one of the phases. Finite-size effects associated with 
the limited size of individual domains, while possibly rel- 
evant, fall outside the scope of the present paper. The 
small radius of the tubes (about 40 nm) does not allow 
a quantitative determination of the shape of the junc- 
tions Pij, and for this reason we restrict ourselves to a 
minimal model which emphasizes the roles of the most 
relevant physical parameters. We show that tube fission 
can be driven either by the line tension or by the jump 
of the elastic coefficients at the interface between the two 
phases, and we compare the two processes that undoubt- 
edly both contribute in the experiments. 

Model. We use an elastic membrane free energy as 
introduced by Canham and Helfrich and numerically 
determine equilibrium junction shapes. Fig. gives a 
schematic representation of the tube and the coordinate 
system used in the following. Our axisymmetric surface is 
parametrised by the arc length s along the contour and 
described by the local tube radius r(s) and the angle 
i/)(s). They are related by the geometric relations r = 
cos?/; (dots denote derivatives with respect to s). The 
interface is located at z = s = 0. 

The free energy of the system is obtained by extending 
the description of tubular membranes [Tjj to the specific 
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FIG. 2: Schematic representation of the junction. 



case of a biphasic tube |2(J, [2l| as follows: 
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The two phases are denoted by a and /?, and for each 
phase i the free energy is integrated over its membrane 
area £V The Ki and k>^q are the bending- and Gaus- 
sian rigidities of the respective phases. This free energy 
includes the bending energy to lowest order in the princi- 
pal curvatures, where H is the mean curvature and K the 
Gaussian curvature. The two layers of the membrane are 
assumed to be symmetric - both phases contain choles- 
terol molecules which have a high flip-flop rate. Any 
stress due to area differences between the leaflets or to 
an asymmetry of the layers is thus quickly relaxed. Fi- 
nally, Lagrange multipliers o~i are introduced to ensure a 
constant surface in each phase. These G{ are interpreted 
as surface tensions. We take our tube to be infinite, and 
assume the presence of a lipid reservoir. In the experi- 
ments, such a reservoir is provided by the large mother 
vesicles from which the tubes are drawn. Provided the 
area per lipid remains constant during the process this 
implies a constant surface tension in each of the phases. 

The interface between the two phases is described by a 
jump in the values of the bending rigidities k!q and in 
the surface tension G{ , and by a positive line tension r at 
the interface dVt. The last term in the free energy is the 
work performed by the external force / needed to pull 
the tube. We neglect the small effect of pressure [lij . 

The variational derivation of the shape equations of 
the surface has been detailed elsewhere |23j, and yields 
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Far away from the junction, we recover homogeneous 
cylindrical tubes with ip = ir/2 and Ri = (^/^c^) 1 / 2 . 
Mechanical equilibrium implies that the forces at both 
extremities are equal and that / = 27r(2cr^) 1 / 2 , which 
imposes that cr a /aj3 = Kpj The surface tension jumps 
discontinuously across the interface. 
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FIG. 3: Numerical shapes of the junction for various line tensions and differences of elastic rigidities in dimensionless units. The 
length scale is the radius of phase a (R a = 1), the energy scale is the bending rigidity of phase a so that K a = 1. (a): shapes 
for various ratios of bending rigidities. The line tension vanishes (r = 0) and the Gaussian rigidities are equal (A^g = 0). The 
values of Kp/hia are 1.25, 1.5, 1.75, and 2.0 (b): shapes for various line tensions. The elastic rigidities are equal: Akg = and 
k& — Kp- The values of the line tension are r = 0.5, 1.0, 1.5 and 2.0. (c): shapes for various differences in Gaussian rigidity. 
The line tension vanishes (r = 0) and the bending rigidities are equal (« a = Kp). The values of the difference in Gaussian 
rigidity are Ak g = n% - = -1-0, 1.0, 2.0 and 4.0. 



The mismatch between constants such as the bend- 
ing rigidities appears only in the boundary conditions 
and strongly affects the interface shape. At the interface 
(s = 0), four boundary conditions must be satisfied. Two 
conditions are the continuity of the radius r(s) and the 
angle ijj(s) poj : two additional conditions stem from the 
variational procedure and relate the first and the second 
derivatives of the angle i/j on each side of the interface to 
the values of r, ^, K a , up, Akq = ^g — anc ^ r - 

Results. Fig. |3| illustrates the different effects that line 
tension and differences in elastic rigidities individually 
have on the two-phase tube. The first possible discon- 
tinuity at the junction is a jump in bending rigidities 
(Fig. ^(a)). The ratio of the bending rigidities in the 
two phases n = Kp/fra fixes the ratio of the radii away 
from the junction and of the surface tensions in the two 
phases. In the absence of both line tension and jump 
in Gaussian rigidity, the radius decreases smoothly from 
the values of the more rigid phase to the value in the less 
rigid phase, but with a remarkable structural feature - a 
small plateau (i.e. a membrane region with a horizontal 
tangent) occurs around the junction. Thisplateau is also 
given by an analytical linear calculation |24j. 

When line tension dominates (Fig. the radius 

at the interface decreases with increasing line tension. 
It vanishes for a huge line tension. Note that our de- 
scription breaks down at scales comparable to the bilayer 
width. Despite the fact that the radius goes to zero the 
mean curvature remains finite; in the highly pinched limit 
a saddle point develops at the neck which keeps the total 
curvature energy finite. 

When the discontinuity in Gaussian rigidities domi- 
nates (Fig. Elfc/), numerical evidence suggests that the 
neck radius does not decrease all the way down to 
zero. Moreover, stability arguments given below impose 
a bound on the maximum absolute value of Akg- How- 
ever, the presence of the neck favors the breaking process. 
In this case, fission does not occur exactly at the interface 



but at the neck. One thus expects to find, after fission, a 
small patch of one phase still attached to the other phase. 
Since details at the length scale of the neck itself cannot 
be resolved experimentally, this effect might be relevant 
to determine the dominant fission mechanism. 

Discussion. For general experimental conditions, all 
three effects are superimposed at the junction. A quanti- 
tative analysis of the shape in order to extract the various 
parameters is then difficult, especially as little to noth- 
ing is experimentally known about the precise shape of 
the junction. Typical values of the bending rigidity of 
liquid bilayers are around 25 ksT, and the rigidity of the 
liquid ordered phases can be up to several times higher. 
Recently, the bending modulus of a heterogeneous vesi- 
cle has been obtained by comparing the experimental 
shape to numerical solutions of the shape equations 25]. 
The Gaussian rigidity Kq* is notoriously difficult to mea- 
sure experimentally, but a recent study cites values of 
Kq = — 0.83/^ j2£|. Stability arguments impose that 
-2Ki < k$ < 0. 

The equilibrium free energy of the tube can be calcu- 
lated from Eq.JIJ and allows a discussion of the stability 
of the tube and of its fission. We show in Fig. 01 the free 
energy of a tube as a function of the dimensionless ra- 
dius at the neck r nec k/^o in the specific case of K a = Kp 
and ^. This energy is maximal for a vanishing 

radius: at this point, the membrane is maximally bent. 
Fission of the tube by pinching requires one to cross this 
energy barrier. The free energy of the ruptured tube is 
also shown on the figure. It is lower than the top of 
the barrier by the contribution of the Gaussian curva- 
ture due to the change in topology upon rupture, which 
equals Attkq. Notice that the bending energy does not 
change upon rupture: at vanishing radii, the neck is a 
saddle point with vanishing mean curvature 0. The 
ruptured tube corresponds to a transient shape since in 
the absence of an applied force, the tubes retract to form 
two spheres. In the absence of line tension the tube is 
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tal precision is not yet at a level where these results can 
be compared directly to our calculated tube shapes, we 
have also considered the breaking time of a two-phase 
tube. Our modeling, based on an energetic approach, 
predicts a strong dependence of the fission dynamics on 
the elastic properties of the phases and yields results that 
are in good agreement with the experimental data. 
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FIG. 4: Schematic energetics of fission. The solid curve plots 
the free energy of a tube pinched by line tension as a function 
of the dimensionless neck radius. Every equilibrium radius r* 
has a corresponding energy = £bend+£r which defines in 

turn an energy barrier for fission £ gap and a free energy gain 
upon fission AJF. For clarity, this figure assumes identical 
elastic rigidities on both sides. 

uniform r(z) = Rq and its energy is zero, and a homoge- 
nous tube is thus thermodynamically stable only if the 
free energy of the ruptured tube is positive. Numerically, 
we have determined this stability limit as k,q > —1.29k. 

The values of the parameters then fix the value of 
^neck/^o- We have also evaluated the energy barrier 
against fission by pinching from this macroscopic model. 
Note, however, that this is only a lower bound to the 
real energy barrier, as it ignores effects at the molec- 
ular lengthscale which certainly is attained when the 
neck becomes very thin. To compare our results to 
the experiments, we have computed the various ener- 
gies at the following (measured or realistic) parame- 
ter values. With bending rigidities n a = 40/c#T = 

1.6 • IO" 19 J, Kp = 70 fc#T = 2.9 • IO" 19 J, Gaussian rigidi- 
ties «g = -33.2 k B T = -1.38 • IO" 19 J, k% = -58 k B T = 
— 2.3-10 -19 J, surface tensions 1-10 -6 N/m in phase a and 

5.7 • 10" 7 N/m in phase /?, and a line tension 7 • 10" 12 N, 
we have determined the height of the energy barrier to 
be £ gap = 7.8 k B T. If we assume that fission is a ther- 
mally activated process [2?|, the average time until fis- 
sion tb occurs is tb = to exp £ gap /&B^- Using a hydro- 
dynamic argument, we estimate the basic time scale as 
to = rjRa/Ka, where n is the viscosity of water. For 
the parameter values cited above this yields a timescale 
to ~ 1.44 • 10 -4 s. We thus expect the experimental time 
until fission to be approximately 350 ms. This is in good 
agreement with the experimentally observed typical time 
for fission, which is of order 1 s. 

Conclusion. We have studied the behavior of a multi- 
phase membrane tube using thermodynamic arguments. 
The shape of the junction between two domains depends 
on three quantities: the line tension of the interface and 
the jumps in the two elastic constants. While experimen- 
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